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ABSTRACT 

Stellar populations in spiral bulges are investigated using the Lick system of spectral 
indices. Long-slit spectroscopic observations of line-strengths and kinematics made 
along the minor axes of four spiral bulges are reported. Comparisons are made be- 
tween central line-strengths in spiral bulges and those in other morphological types 
(elliptical, spheroidal (Sph) and SO). The bulges investigated are found to have central 
line-strengths comparable with those of single stellar populations of approximately so- 
lar abundance or above. Negative radial gradients are observed in line-strengths, sim- 
ilar to those exhibited by elliptical galaxies. The bulge data are also consistent with 
correlations between Mg 2 , Mg 2 gradient and central velocity dispersion observed in 
elliptical galaxies. In contrast to elliptical galaxies, central line-stren gths lie within the 
loci defining the range of <Fe> and Mg2 achieved by Worthey's ( 1994 ) solar abun- 
dance ratio, single stellar populations (SSPs). The implication of solar abundance 
ratios indicates significant differences in the star formation histories of spiral bulges 
and elliptical galaxies. A "si ngle zone with in-fall" model of galactic chemical evolu- 
tion, using Worthey's (1994) SSPs, is used to constrain the possible star formation 
histories of our sample. We show that the <Fe>, Mg2 and H/3 line-strengths observed 
in these bulges cannot be reproduced using primordial collapse models of formation 
but can be reproduced by models with extended in-fall of gas and star formation (2- 
17 Gyr) in the region modelled. One galaxy (NGC 5689) shows a central population 
with a luminosity weighted average age of ~ 5 Gyr, supporting the idea of extended 
star formation. Kinematic substructure, possibly associated with a central spike in 
metallicity, is observed at the centre of the Sa galaxy NGC 3623. 

Key words: galaxies: abundances, galaxies: evolution, galaxies: formation, galaxies: 
spiral, galaxies: stellar content. 



1 INTRODUCTION 

The recently dev e loped Lick system of spectral indices 
(Faber et al. 1985; Worthey et al. 1994) has provided pow- 
erful tools for the investigation of composite populations. 
Diagnostic plots, such as <Fe> vs H/3 and Mg2 vs <Fe> , 
(where <Fe> is the average of Fe5270 and Fe5335 indices) 
break the age /metallicity degeneracy and have illustrated 
characteristics and trends in elliptical galaxies that constrain 
their possible star formation histories (SFH). This paper 
compares the line-strengths in four spiral bulges with those 
observed in other morphologies and places constraints on 
possible SFHs in the sample. Comparison of the SFH of spi- 
ral bulges to that of elliptical galaxies is informative as the 
two morphologies show numerous similarities. Bulges exhibit 
similar surface brightness profiles and kinematic structures 
to ellipticals and form an overlapping co ntinuation of the el- 
liptical locus in the Fundamental Plane (Bender, Burstcin & 



ever, despite similarities in physical and photometric prop- 
erties, the two galaxy types are found in very different envi- 
ronments, with the majority of ellipticals located in clusters, 
while spiral galaxies are preferentially found in the field. An- 
other important factor in the environment of spiral bulges 
is the presence of the disc, which severely limits possible 
merger histories. 



Observations of elliptical galaxies have shown them to 
possess high central line-strengths. Gradients are observed in 
metallicity sensitive features indicative of dissipation in the 
formation process. The H/3 index shows no indication of sys- 
tematic gradients in age (Carollo, Danziger & Buson 1993; 
Davies, Sadler & Peletier 1993; Fisher, Franx & Illingworth 
1996; Gorgas et al. 1997; Vazdekis et al. 1997). Modelling of 
Lick indices in elliptical galaxies, for comparison t o observa- 



tion, has been carried out by a numbe r of authors ( Worthey. 



Faber 1 99£ ) . Colours and colou r gradients hi t he two types Dfflmia^k_Jojies_^996 Srcggio 1997 



Vazdekis et al. 1997 



of spheroid are also similar ( Balcells fc Peletier 1994 ). How- Sansom 1- Proctor 199^ ). These studies have shown that the 
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high line-strengths found at the centre of elliptical galaxies 
can only be achieved if the bulk of the observed population 
forms from pre-enriched gas. Two possible processes have 
been proposed to produce this pre-enrichment. One possi- 
bility is an initial mass function (IMF) biased towards high 
mass stars in the early stages of galaxy formation. In this 
way the first few generations of stars generate large quan- 



titie s of metals (Gibson & Matteucci 1997; Vazdekis et al 
199 7|). j \lternatively, a period of star formation (SF), and 
thus interstellar medium (ISM) enrichment, with a Salpeter 
IMF prior to the main burst, can produce the necessary 
pre-enrichment. This delayed burst of SF may correspond 
t o a merger /interaction m odel of elliptical galaxy formation 
( |3ansom fc Proctor 1998b - 

The second important feature of line-strengths in el- 
lipticals is their high central Mg2 with respect to <Fe> 
when compared to models assuming solar abundance r atios 
( IWorthcy, F a bcr fc Gonzalez 1992|; Pavies et al. 199$ |Gor 



gas et a |. 1997] ; Greggio 1997). This Mg2 excess is interpreted 



as indicating [Mg/Fe] greater than solar (Davies et al. 1993) 



High [Mg/Fe] can be achieved by star formation times of & 
1 Gyr, in accordance with the s hort time-scales of gas in - 
flow in both merger/interaction (Barnes & Hernquist 1996) 



and primordial collapse models of elliptical galaxy formation 
(Theis, Burkert & Hensler, 1992). Short bursts of star forma- 
tion ensure that the metal enrichment of the ISM from which 
the bulk of the population forms, is mainly produced in the 
Type II supernovae that are the main sites for production 
of Mg. The biased IMF proposal also naturally produces an 
excess of Type II supernovae (and consequently Mg2) due to 
an over-production of high mass stars. While both merger 
and primordial collapse (with biased IMF) models can suc- 
cessfully reproduce the indices in the Mg2 vs <Fe> plane, 
the high central H/3 values of many elliptical galaxies are a 
challenge to primordial collapse models. In contrast, merger 
models clearly predict the existence of younger populations 
in these objects. 

Observation s of other galaxy mo rphologies, s uch as SOs 
flSU'chenko 1993) ; [Fisher et al. 1996| ) and Sphs ( |Gorgas et 
al. 199 <|), show differences between central line-strengths in 
various galaxy types in the <Fe> vs Mg2 plane. SOs tend 
to exhibit weaker Mg2 for a given <Fc> than do Es while 
still showing an excess with respec t to v alues for solar abun- 
dance ratio SSPs. Gorgas et al. ( 1997 ) showed that Sphs 
roughly follow solar abundance ratio SSP loci in the <Fe> 
vs Mg2 diagnostic plot. Low luminosit y spi ral bulges ob- 
served by Jablonka, Martin & Arimoto (199£) also span the 
region of the Mg2 vs <Fe > plot occupied by Worthey's so- 
lar abundance ratio SSPs ( Worthey 199§| ), while some high 
luminosity bulges in their sample exhibit greater than solar 
[Mg/Fe]. 

Observations of individual stars in the bulge of our own 
Galaxy, while hampered by obscuration, have provided in- 
teresting insights into the population present. Minniti et al. 
(1992) report a broad range of metallicities (—2 <[Fe/H]< 



+ 1) in K giants in the Galactic bulge and indicate the pres- 
ence of a metallicity gradient. Observations of resolved M 
giants in the centre of the bulge of our own Galaxy indicate 
approximately solar abunda nces and also indicate a gradient 
in metallicity (Frogel 1998). Recent episodes of star forma- 
tion at the centre of the bulge, with several bursts having 
occurred over the past few 100 Myr, have also been sug- 



gested (Blum, Sellgren & Depoy 1996; Frogel 1998). These 
observations suggest that the population of the bulge of our 
own Galaxy has a complex and extended star formation his- 
tory. 

Our long-slit observations of spiral bulges allow compar- 
ison of central line-strengths and gradients to those in other 
galaxy types and constrain possible SFH in the bulges. We 
probe whether the similarities between bulges and elliptical 
galaxies in morphological and photometric properties are the 
result of similar formation processes, or are rather the result 
of the evolution of differing histories to similar present-day 
forms. In Section ^ the observations, data reduction and er- 
ror estimates are outlined. Results of our analysis and com- 
parison with other Hubble types are reported in section ^[ 
Comparison to models are detailed in section Our conclu- 
sions and the direction of future work are outlined in section 

1 



2 OBSERVATIONS & DATA REDUCTIONS. 
2.1 Observations. 

Using the Palomar 5.08 meter telescope on 1995 March 31, 
long-slit spectroscopy was carried out along the minor axes 
of the bulges of 4, approximately edge-on, spiral galaxies 
(table and fig. [j]). We chose to observe edge-on systems 
along their minor axes as disc obscuration of bulge light is 
thus concentrated on one side of the bulge. A fifth galaxy 
(NGC 4013) was also observed but appeared to have a bright 
foreground star slightly offset from the centre of the galaxy. 
Consequently, this galaxy was excluded fro m data reduc- 
tions. The Double Spectrograph was used (Oke fc Gunn 



1982) with a 5700 A dichroic filter and Texas Instruments 



CCDs on both red and blue arms. Wavelength coverage in 
the blue was 4560 - 5430 A at a resolution of 1.1 A per pixel. 
In the near infra-red the range 8235 - 8860 A was observed 
at a resolution of 0.8 A per pixel. The blue wavelength range 
was chosen to include the H/3, M g and Fe absorption features 
calibrated by Fabe r et al. (|l985 [ ) and extensively observed in 
elliptical galaxies (|Davies et al. 1993; Hes & Peletier 1993|; 



Fisher et al. 1996 



|Gorgas et al. 1997 



Vazdekis et al. 1997) 



The red wavelength range spans the Call IR triplet and was 
intended for the study of both indices and kinematics. Spa- 
tial resolutions of the CCDs were 0.78 and 0.58 arcsec per 
pixel on blue and red arm respectively. Seeing was ~ 1.5 
arcsec and a 1 arcsec slit was used. Flux and velocity stan- 
dard stars were also observed. Object frames were alternated 
with arc lamp exposures to allow for accurate wavelength 
calibrations. Frames of a continuum source (tungsten lamp) 
were also obtained for the purpose of flat-fielding. Galaxy 
exposure times of 5900, 7200, 4200 and 1200 seconds were 
obtained for NGC 2654, NGC 3623, NGC 4565 and NGC 
5689 respectively. A maximum single exposure time of 2400 
seconds was adopted to facilitate cosmic ray removal. 



2.2 Data Reduction. 

Unless otherwise stated, data reduction was carried out us- 
ing the FIGARO package of Starlink software. A bias level 
from the over-scan region was subtracted from each frame. 
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Table 1. Data for four spiral bulges observed with the Palomar 5m telescope. Hubble Type and B magnitudes are from de Vaucoulcurs 
et al. (1991) (hereafter RC3). Inclinations are from Guthrie (1992). Distances are from Tully (1988) assuming H =75 kms _1 Mpc _1 . 
Recession velocities given are from RC3 (HI data) and this work (optical data). The velocity dispersions given (r ry) are the averages 
across all data points. This is the value used for calculation of broadening required to the Lick system (see section 2.3). 



Galaxy 


Hubble 


Inclination 


Distance 


B 


Recession Velocity 


cry 


Scale 




Type 


(deg) 


(Mpc) 


(mag) 


RC3/This work 


(kms ) 


(pc/arcsec) 












(kms" 1 ) 






NGC 2654 


Sab 


90 


23.3 


12.8 


1341/1320 


136 


110 


NGC 3623 


Sa 


79 


7.3 


9.6 


807/809 


152 


35 


NGC 4565 


Sb 


90 


9.7 


10.3 


1227/1235 


155 


16 


NGC 5689 


SO/a 


78 


35.6 


12.7 


2160/2221 


149 


168 



N 



E 
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Figure 1. Digital Sky Survey images of galaxies with slit positions along minor axes. The slit length is 2 arcmin in all images. Slit 
position was determined by assuming the slit passes through the telescope pointing position along the minor axis. The location of the 
bulge within the slit was determined by ensuring that the peak of the luminosity profile on the slit matched the position of the maximum 
count rate on the CCD. Upper left: NGC 2654. Note that the slit appears off-centre by approximately 5 arcsec. Upper right: NGC 
3623. Sky is not reached on either side of the galaxy. However, sky was estimated from extreme edges of slit (see section Bjh . Lower 
left: NGC 4565. Sky is not reached on the side of the galaxy most affected by the disc. However, a reasonable sky estimate is achieved 
on the side least affected by disc. Lower right: NGC 5689. Note the slit appears off-centre by approximately 5 arcsec. 



Due to the presence of localised features in the flat-field 
frames, each object was divided by the normalised, average 
flat-field. This biases the data with the spectral response of 
the tungsten lamp. However, this is smooth and is removed 
during flux calibration. Cosmic rays were removed by in- 
terpolation between adjacent pixels. Wavelength calibration 
was achieved by interpolating between arc lamp exposures 



bracketing galaxy frames and was accurate to ~ 0.1 A. Spec- 
tra were then flux calibrated and extinction corrected. For 
NGC 2654, NGC 3623 and NGC 5689 sky subtraction was 
carried out by interpolation across the galaxy. For NGC 3623 
sky was not fully attained on either side of the galaxy. This 
results in a small contribution of disc light to the sky es- 
timate. The error associated with this sky was estimated 
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from the 1 sigma variation in index values when sky esti- 
mates from successive galaxy exposures (NGC 2654, NGC 
3623 and NGC 4565) were used. For NGC 4565 sky was 
only reached along the slit on one side of the bulge. In this 
case the sky was estimated from just this side of the galaxy. 
Variati on of the sky across the slit was shown to be less than 



10% of iky level for frames clearly reaching sky on both sides. 



Finally, galaxy frames were shifted and co-added to form a 
single frame for each galaxy. 



2.3 Measurement of Kinematics and Indices. 

Kinematics were measured in order to correct our data to 
the dispersion of the Lick system. Measurements of kinemat- 
ics were carried out on both red and blue data using both 
cross-correlation and Fourier quotient techniques within the 
IRAF software package. Results of the two techniques were 
in agreement to within 5 kms -1 . However, the Fourier quo- 
tient technique was used as the statistical errors associated 
with cross-correlation were substantially larger due to the 
insensitivity of this method at the low velocity dispersions 
typical of spiral bulges. As velocity dispersions in our sam- 
ple are all less than the dispersion required for calibration to 
the Lick system, the additional broadening required ( as) 
can be calculated (assuming Gaussian profiles) by; 



star is not suitable for estimation of correction factors to the 
Lick system. Lack of sufficient data to fully calibrate to the 
Lick system represents a limitation on the Palomar data set. 
However, corrections to the narrow line features discussed 
here (Fe4 668, H/3, Fe5270 and Fe5335) have been shown to 
be small (|Gorgas et al. 1997 ; fVazdekis et al. 1997| ; [Worthey 
; Ottaviani 1997). For Jones ' ( |1996| ) observations correction 



„2 2 „2 

a L - a j - a v 



The Lick resolution (<tl) is approximately constant at 3.57A 
over the 4500-5500 A rang e covered by our observations 
(Worthey & Ottaviani 1997). Instrumental broadening (crj) 
was estimated from arc lines and found to be 1.12 A. Veloc- 
ity dispersion broadening of the galaxy (av) was estimated 
as the average over the whole galaxy in the blue (table [j]). 

After broadening the spectra to the Lick resolution, the 
indices were evaluated by our own code, using the band 
wavelength range definitions supplied by G. Worthey on his 
home page. Worthey also kindly provides data with which 
our code could be tested. Differences between the Worthey 
derived indices and our own for the provided spectra were 
< 0.03 A for the line features and < 0.002 mag for the Mgi 
and Mg2 molecular bands. These discrepancies are smaller 
than differences caused by re-calibration of Worthey's data 
to our wavelength resolution and are probably the result of 
the effects of differences in the handling of partial bins. They 
are also significantly smaller than the total errors (section 
2.4). A more important source of uncertainty in our data is 
the conversion to the Lick system (which was not flux cali- 
brated). One of the observed velocity standards, HD 139669 
is also one of the Lick calibration stars. Indices derived from 



our observations of this star and the data of Jones (1996) 



(data available on AAS C D-ROM Series Vol. VII) were com- 
pared to the Gorgas et al. (1993) published data. Differences 
between indices derived from our data and those of Jones 
are typically within errors. However, there is a large dis- 
crepancy between our measurement and the published Lick 
data for the H/3 and Fe5015 indices (also reproduced in the 
data of Jones) which is probably caused by unpredictable 
wavelength shifts in the Lick IDS spectra (private commu- 
nication with J. Gorgas). Differences between our measured 
indices and the published Lick data for this star are given 
in table |[ HD 139669 (HR 5826) al so lie s well away from 
the mean in fig. 10 of Gorgas et al. (1993). Therefore, this 



to the flux calibration sensitive Mg2 index was found to be 
small (0.01 mag increase in the measured Mg2 index). The 
good agreement between our measurement of the Mg2 index 
of HD 139669 and that of Jones suggests that our correction 
factor to the Lick system should be of similar magnitude (~ 
+0.01 mag), although we have not applied this correction. 

Initial measurements of indices showed unusual gradi- 
ent inversions in the Call, Mgi and Mg2 indices at the cen- 
tres of the bulges (defined as the peak in brightness) . These 
were found to be the result of focus variation caused by 
distortion of the CCD surface. This effect redistributes flux 
spatially between bins, introducing false continuum varia- 
tions. Call and molecular band indices are sensitive to this 
type of distortion due to the large separation of their side- 
bands. The problem was most pronounced in the red (this 
CCD has subsequently been replaced at Palomar). As a re- 
sult, line-strength measurements were not carried out on 
the red arm data. Kinematic analysis were still possible as 
the continuum shape is removed prior to line-broadening 
measurement. Stellar data from the blue arm had a smooth 
point source FWHM variation of <20% across the wave- 
length range. Even so, features are still introduced in the 
Mgi and Mg2 indices in the centre of the galaxies where 
luminosity gradients are greatest. A minimum binning of ~ 
1.5 arcsec has minimised this effect and matched the see- 
ing conditions. Central values for purposes of modelling and 
comparison to ellipticals were calculated using the central 3 
arcsec (see section [j.2] ). 

2.4 Estimation of Errors. 

Detailed analysis of errors affecting NGC 3623 are shown 
in table ^| for illustration. For all galaxies, binning of data 
was selected to maintain an approximately constant Poisson 
error outside the central regions. For NGC 2654, NGC 4565 
and NGC 5689, sky errors were evaluated by varying the 
region of sky used for sky estimation. The 1 a variation in 
estimated index value is taken as the error. For NGC 3623, 
which failed to reach sky on both sides of the galaxy, sky 
error was estimated as the 1 sigma variation in indices de- 
rived utilising sky estimates from NGC 3623 and the brack- 
eting observations of NGC 2654 and NGC 4565. Variation 
of both velocity dispersion and recession velocity about the 
assumed average values (table Q) were both ~ ± 35 kms" 1 . 
These variations are significantly larger than errors associ- 
ated with wavelength calibration and statistical errors from 
the Fourier quotient technique (~ 9 kms" 1 ). Consequently, 
velocity dispersion errors were estimated as the variation in 
indices caused by increasing the galaxy velocity dispersion 
assumed in calculation of broadening to the Lick system by 
35 kms" 1 . Index errors due to recession velocity uncertainty 
were estimated as the variation in indices caused by increas- 
ing the red-shift estimate by 35 kms" 1 . 

The magnitude of the sky error varied between galax- 
ies. Consequently, unless otherwise stated, figures and tables 
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Table 2. Error table for NGC 3623 Above: Index values for the central 3 arcsec and the associated Poisson noise are given. Sky errors 



are derived as detailed in section 



2.4. Calculations of both velocity dispersion and recession velocity errors (a and V respectively) assume 
l 



approximate velocity uncertainties of ± 35 kms . Total errors are Poisson, sky, a and V errors added in quadrature. Comparison to 
HD 139669 represents published Lick values minus values obtained from our observation of this star. These offsets are not included in 
the calculation of total error due to their uncertainty. Below: Error table for outer regions (~ 8 arcsec from centre) of NGC 3623. 



CENTRAL VALUES 



Index 


Index 


Poisson 


Sky 


a 


V 


Total 


Comparison 




Value 


Noise 


Error 


Error 


Error 


Error 


to HD 139669 


Fe4668 


8.64 A 


0.07 


0.05 


0.05 


0.06 


0.12 


0.15 


H/3 


1.29 A 


0.03 


<0.01 


0.03 


0.02 


0.05 


0.60 


Fe5015 


4.91 A 


0.07 


0.07 


0.23 


0.04 


0.25 


-0.93 


Mgb 


4.51 A 


0.04 


0.02 


0.09 


0.02 


0.10 


-0.03 


Fe5270 


3.31 A 


0.04 


0.01 


0.10 


0.01 


0.11 


-0.47 


Fe5335 


3.07 A 


0.05 


<0.01 


0.17 


0.02 


0.18 


0.14 


Mgi 


0.133 mag 


0.002 


<0.001 


0.001 


< 0.001 


0.002 


0.014 


Mg 2 


0.279 mag 


0.002 


<0.001 


0.001 


0.001 


0.002 


0.016 


OUTER VALUES 


Fe4668 


6.88 A 


0.12 


0.28 


0.07 


0.05 


0.32 


0.15 


H/3 


1.22 A 


0.06 


0.01 


0.01 


0.01 


0.06 


0.60 


Fe5015 


3.47 A 


0.12 


0.47 


0.21 


0.11 


0.54 


-0.93 


Mgb 


4.12 A 


0.07 


0.12 


0.07 


0.01 


0.16 


-0.03 


Fe5270 


2.79 A 


0.08 


0.03 


0.09 


0.01 


0.12 


-0.47 


Fe5335 


2.59 A 


0.09 


0.02 


0.13 


0.01 


0.16 


0.14 


Mgi 


0.107 mag 


0.003 


0.001 


0.001 


< 0.001 


0.003 


0.014 


Mg 2 


0.247 mag 


0.004 


0.003 


0.001 


0.001 


0.005 


0.016 



show combined Poisson and sky errors. Velocity dispersion 
and recession velocity errors for other galaxies were simi- 
lar in magnitude to those in NGC 3623. In the centres of 
the bulges Poisson noise and velocity dispersion errors dom- 
inate. While for most cases sky error dominates in the outer 
regions, total errors in Fe5270 and Fe5335 continue to have 
significant contributions from velocity dispersion error. Gra- 
dient errors given in table ^ were estimated as the statistical 
errors associated with linear fitting. 



3 RESULTS. 
3.1 Kinematics. 

Results of recession velocity and velocity dispersion analy- 
sis are shown in fig. || Galaxy centres are defined as the 
position of the luminosity peak along the slit. Reasonable 
agreement exists between results for red and blue data with 
average values for both recession velocity and velocity dis- 
persion differing by less than 20 kms -1 . However, systematic 
differences are present in the centres of NGC 3623 and NGC 
4565. Central recession velocities derived from the red data 



are higl ly susceptible to systematic error due to focus varia- 
tion (section |2.3[ ). Consequently, analysis of indices has been 
based on the blue recession velocity data only. In NGC 4565 
(a known LINER, Keel, 1983), the maximum in the blue ve- 
locity dispersion profile is offset by ~ 5 arcsec with respect 
to both the luminosity and metallicity peaks (see figs ^ and 
|i|). This would appear not to be an internal extinction effect 
as metallicity continues to rise with luminosity to the centre. 
Photometric observations in the near-infrared of NGC 4565 



in the centre, show a flattening of the colour gradient that is 
evident along the minor axis at distances greater than 5 arc- 
sec (Rice et al. 1996; fig. 2a). The change in colour gradient 
is coincident with the peak in blue velocity dispersion data 
in fig. ^| and and the edge of the broad metallicity peak in 
fig. [4 The central region of radius 5 arcsec is also the emis- 
sion region of NGC 4565. Strong [OIII] 5007 A and 4959 A 
emission lines are present in this region affecting the Fe5015 
index, while the H/3 index exhibits features due to emission 
line filling (fig. ||). Thus, the data suggest that the core of 
NGC 4565 contains a blue, metal rich, kinematically cold 
sub-population perhaps associated with ongoing star forma- 
tion. The presence of such a population would also cause 
a displacement between red and blue luminosity peaks as 
observed. 

A dip in velocity dispersion is seen in the centre of NGC 
3623 and is reproduced in both red and blue data sets. Such 
velocity dispersion profiles have been observed along the ma- 
jor axis of spiral galaxies (e.g. Bertola et al. 1996; Bottema & 
Gerriston 1997) . Kinematic models that reproduce such dips 
in velocity dispersion include the presence of a kinematically 
distinct comp onent within the bulge such as an iso thermal 



nuclear cor e ( Bottema fc Gerritsen 1997 ), a disc (Bertola 
et al. 1996) or a bar flFriedli 1996[ ). Chemo-dynamical mod- 



elling of the e ffect of the presence of a bar on disk galaxies 
(Friedli 1998) has shown that metallicity gradients are re- 



( Rice et al. 1996 ) , rather than indicating reddening by dust 



duced in the plane of the bar around and outside the region 
of the co-rotation radius, while steeper gradients can occur 
at radii well within co-rotation. A flattening of the metallic- 
ity gradients is also induced along the minor axis. The gra- 
dients in metallicity sensitive indices along the minor axis of 
NGC 3623 (figs|), which flatten rapidly away from the cen- 
tre, are consistent with such models. Friedli's ( 1994 ) models 
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Figure 2. Spatially resolved kinematics along the minor axes of 4 spiral bulges. Recession velocities and velocity dispersions from 
both the blue and red Palomar data are shown. Both sets of data were derived by the Fourier quotient technique. Error bars represent 
statistical errors of this technique and do not include systematic errors. 



also predict inflow of gas within the co-rotation radius re- 
sulting in a central burst of metal rich star formation. The 
sharp metallicity peak, evident in both Fe and Mg indices in 
fig [s| may be the result of this process. Photometric obser- 
vations by Burkhead et al. 1981 and recent Hubble archive 
images of NGC 3623 show boxy central isophotes which may 
also be indicative of the presence of a bar. Alternatively, the 
kinematic feature may be the result of an embedded stellar 
disc. If this is the case, either the disc has a scale-height 
of ~ 150 pc (major axis disc) or the disc rotates along the 

minor avis with a radius of at 15Q pc Hnwpwr t.hprp is 



3.2 Lick Indices. 

In this section we present the results of index evaluations 
and plot the data on several diagnostic plots for compar- 
ison to models of SSPs and observations of other galactic 
morphologies. Index measurements across the 4 bulges are 
shown in table "j"] and fig. [j. Negative radial gradients can 
be seen in all the metallicity sensitive indices with the pos- 
sible exception of Fe5015, which is strongly affected by the 
[OIII] 5007 A emission line. Mgb is affected by emis sion of 
the [Nil 5199 A double t in one of the side-bands flGoud 



no evidence of the kinematic signature of a minor axis disc frooi j & Emsellem 1999), However, the enhancement of this 



on this scale in the recession velocity measurements (fig. |2j). 
On balance it would seem likely that the centre of this spiral 
bulge harbours a small scale bar. (Note also that RC3 gives 
the classification of NGC 3623 as SXT1. This indicates a 
"mixed" bar/no bar type). 

As all the indices described here are from the blue data, 
and the red data suffer degradation by focus variation, we 
have used the blue velocity dispersion and recession velocity 
values in our analysis of the indices. Average values across 
the galaxies were used for both recession velocity and veloc- 
ity dispersion in recognition of the systematic uncertainties. 
The resultant recession velocities are in good agreement with 
values given in RC3 (table m) . 



index caused by emission appears to be small in these data 
for all galaxies, with the possible exception of the LINER 
NGC 4565. The age sensitive H/3 index can also suffer from 
emission line filling, particularly near the centres of bulges. 
This is most clearly seen in NGC 4565 (fig. []). Dips in H/3 
at the centre of this galaxy are associated with two separate 
emission regions and are reflected by features in the Fe5015 
index. The dip in H/3 in the central region of NGC 3623 may 
also indicate weak emission. However, this galaxy exhibits 
sharp central peaks in metallicity sensitive features which 
coin cide with a dip in the velocity dispersion (see section 
3.1). The combined effects of the the influence of the metal- 
licity peak and emission on the H/3 index in the centre of 
this galaxy makes interpretation difficult. The symmetry of 
indices about galaxy centres indicate that only in NGC 4565 
is disc obscuration significant. 

Fig. ^ shows central values of Mg2 vs <Fc> for our four 
spiral bulges. Also plotted are central index v alues of ellip- 
tical (E) and Sph galaxies (|Gorgas et al. 19*97 ), SO galaxies 
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Figure 3. Indices across the minor axes of four spiral bulges. Radial gradients can be seen in most metallicity sensitive features. The 
central dips in H/3 in NGC 4565 are the result of emission. Fe5015 is also emission affected. The errors bars shown include Poisson noise 
and sky error only. Bin sizes were selected so as to maintain approximately constant Poisson error values beyond the central pixels. The 
increasing influence of the sky error at large radii can be seen. Note the sharp central peak in metallicity sensitive indices in NGC 3623. 
This peak coincides with a dip in central velocity dispersion. 
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Table 3. T op: Central 3 arcsec index values (except H/3 for NGC 4565 where the average of outer values is used to avoid emission: see 
section |3,2| ). Errors (given in brackets) are Poisson noise and sky errors combined. Middle: Index gradients calculated from all available 
data points. For NGC 4565, errors are statistical errors (associated with linear fitting technique) combined with uncertainty in position of 
the galaxy centre. Positional uncertainties are small in other galaxies, consequently, for these, quoted errors are purely statistical errors. 
Bottom: Central velocity dispersions. Quoted errors are statistical errors from the Fourier quotient technique. 1 Peak value 175 kms -1 
at 5 arcsec off-centre. 2 Peak value 190 kms -1 at 5 arcsec off-centre. 
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Figure 4. Mg2 vs <Fe> for elliptical and Sph galaxies (Gorgas et al. 1997) and SO galaxies (Fisher et al. 1996). Results for the 4 spiral 
bulges presented here and the bulge of the Sa galaxy M104 (Hes & Peletier 1993) are also shown. The dotted lines indicate the extremes 
of Worthey's (1994) SSPs and also the solar metallicity locus (Z=0.0). The solid lines indicate the range of index values achieved by 
primordial collapse models with a range of star formation efficiencies from 0.2 to 4.0 Gyr -1 (described in section Dashed lines indicate 
index values achieved by a range of extended inflow models with star formation efficiencies of 4.0 Gyr -1 over 4 to 6 Gyr. Ages of 5, 8, 
12 and 17 Gyr are shown. All models plotted assume solar abundance ratios. 
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(Fisher ct al. 1996 ) and the bulge of M104 (]Hes fc Peletiei] Bell 1995 ), this may indicate either an abundance ratio dif- 



1993;). jrhis diagram illustrates trends in the relative abun 
dance and abundance ratio of different morphological types. 
Dotted lines indicate the extent of the space occupied by 
Worthey's solar abundance ratio SSPs and the solar metal- 
licity locus. The high central strengths of the <Fe> index 
found in E and SO galaxies are also apparent in spiral bulges, 
with bulge indices comparable to those of SSPs with so- 
lar metallicity or above. Fig. ^ also demonstrates that the 
cores of Es exhibit Mg2 values in excess of those possible for 
SSPs with Salpeter IMF and solar abundance ratios. This 
is the well known Mg excess in elliptical galaxies discussed 
in section |l]. SO galaxies lie intermediate between Es and 
the SSPs. The position of the four spiral bulges is consistent 
with solar abundanc e rati os. This is in agre emen t with the 
result of Sil'chenko ( |l993| ) an d Wo rthey's ( [l998| ) interpre- 
tation of the Jablonka et al. ( 1996| ) data for bulges. Outer 
regions of the bulges are also consistent with solar abun- 
dance ratio SSPs (Sansom, Proctor fc Reid 1998). Analysis 



of t he fl ux calibrated data of Jones (Worthey & Ottaviani 



1997) indicated additive corrections of +0.01 mag and -0.06 
A for calibration of Mg2 and <Fe> features respectively to 
the Lick system. While reasonable values for this correction 
for the Palomar data may move some of the spiral bulge 
data points just outside the region of solar ratio SSPs, the 
conclusion that Mg is significantly under-abundant in spiral 
bulges, when compa red to ellipticals, rem ains unchanged. 
The bulge of M104 ( [lies fc Peletier 1993| ) does not appear 
to follow the trend suggested by our sample (fig. ^]), but in- 
stead falls amongst the elliptical galaxies. Therefore, these 
data support the suggestion by Hes and Peletier, that the 
bulge of M104 (classified as SO/a in RC3) more closely re- 
sembles an elliptical galaxy than a normal spiral bulge. 



In fig. g we plot <Fe> against H/3. This diagram pro- 
vides evidence relating to the relative ages of stellar popu- 
lations. The correction factors to transform flux calibrated 
H/3 data to the Lick system have been shown to be small 
(Worthey & Ottaviani 1997). Index values from our data 
are from the central 3 arcsec, except for NGC 4565 where 
the average outer H/3 value was used (omitting central values 
where this index is emission affected). We therefore assume 
no significant gradient in H/3 for NGC 4565. This is consis- 



tent with studies of elliptical galaxies (garollo et al. 1993 



Davies et al. 199S 



Gorgas et al. 1997 



Fisher et al. 1995; Fisher et al. 1996 



Vazdekis et al. 19971) that show no de- 



tectable age gradients. The possible presence of emission in 
the bulges creates uncertainty in the reliability of our H/3 
determination, making interpretation of this index difficult. 
However, it should be noted that NGC 5689 (with H/3 ~ 
2.0 A) shows no sign of central emission which, in any case, 
would only increase the estimate of H/3 for the underlying 
population. Thus, the strong H/3 in this galaxy is compelling 
evidence of an intermediate age (< 5 Gyr) population. 



Fig. ^| locates the spiral bulges in the Fe4668 vs H/3 
plane. The Fe4 668 index is extremely metallicity sensitive 
(Worthey 1994). The correction factor to transf orm Fe4668 
to the Lick syst em has been shown to be small ( |Worthey & 
Ottavia ni 1997 ). Comparison of figs ^ and ^ illustrates that 
higher metallicities are suggested by Fe4668 tha n by <Fe>. 
As the Fe4668 index is highly carbon sensitive (Tripicco & 



ference or a problem with the calibration of Fe4668. The 
main source of carbon enrichment of the ISM is still contro- 



versial. Gustafsson et al. (1999) report falling [C/Fe] with 
increasing [Fe/H] in the solar neighbourhood. They suggest 
that high mass stars are the major contributors of C to the 
ISM. However, both intermediate and low mass stars are 
also known to produce C. Consequently, the interpretation 
of a C over-abundance is difficult. It must also be noted that 
at high metallicities the value of this index changes rapidly 
with metallicity. It is possible that the seeming Fe4668 ex- 
cess is in fact the result of the lack of calibrated SSP data at 
metallicities greater than [Fe/H] = 0.5. Caution is advised 
in interpreting this index at high metallicities. 

3.3 Gradients in Indices and Correlations. 

Gradients were estimated by least squares fitting for the 
<Fe>, Fe4668, Mg2 and H/3 indices plotted versus log radius 
from the galaxy centre (fig. Q) . Data from both sides of the 
galaxies were included in these fits. The gradients obtained 
are given in table |S| Central galaxy regions are marginally 
affected by seeing. However, central index values were in- 
cluded in the gradient estimates due to the limited number 
of data points. As gradients are unaffected by a constant 
off-set to convert to the Lick system, a direct comparison 
with elliptical galaxies can be made. Gradients, central Mg2 
values and central velocity dispersions (oo; table ^) of the 
bulges are consi stent with the correlations reported by Car- 
ollo et al. (1993) in elliptical galaxies. The bulges show Mg2 
gradients of similar magnitude to those in elliptical galaxies 
possessing the same central velocity dispersions and central 
Mg2 values. The gradients are also simila r in m agnitude to 
those reported in Sansom, Peace & Dodd (1994) for 2 bulges 



(NGC 3190 and NGC 1023). The central velocity disper- 
sions and the central index values of our Palomar sample, 
indicating solar abundance ratios, are a lso consistent with 
the pattern suggested by Worthey (1998) that all spheroids 
with velocity dispersions less than 225 kms -1 possess solar 
[Mg/Fe]. 



4 MODELLING CENTRAL 
LINE-STRENGTHS. 

In order to investigate possible SFHs, we use our galac- 
tic chem ical e volution (GCE) code detailed in Sansom & 
Proctor (1998). Briefly, the model calculates star formation 



rate (SFR) and the metallicity of the gas at each time- 
step throughout the history of the region being modelled. 
As all stars formed in each time step are of the same age 
and metallicity, they constitute an SSP. Composite indices 
are then calculated as the luminosity weighted sum of in- 
dices from inter polations between tabulated SSP values of 
Worthey (1994) which cover populations of age 1.5 to 17 



Gyr. The code currently models a single zone allowing for 
gas in-fall. The rate of gas inflow can be varied and can be 
either primordial or enriched to the same level as the gas al- 
ready in the region (simulating inflow from a neighbouring 
region with similar SFH). SFR is assumed to be propor- 
tional to some power (a) of the gas mass density (p) ; SFR 
= Cp a , where C is the star formation efficiency (Gyr -1 ) if 
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a=l. Kennicutt (1989) shows that the value of a in galactic collapse (Carlberg 1984) show that the period of gas inflow 



discs lies between 1 and 2. In all models presented here we 
assume a=l and inflow is enriched. The differences in in- 
dices caused by assuming a = 2 were shown to be small and 
make no significant difference to the conclusions for the mod- 
els presented here. The code permits the modification of the 
star formation efficiency and inflow rates at two points in the 
history of the galaxy. In this way various formation scenar- 
ios can be modelled. The star formation efficiency modelled 
ranges between values for typical Sb galaxies (C = 0.2) up to 
those for star-burst galaxy types (C = 4) . In order to follow 
the enrichment of the gas by feedback mechanisms, such as 
SNII, SNIa and mass-loss from intermediate mass stars, the 
total metal content of the gas is traced in a self consistent 
manner. Data for S NII ( high mass st ars) w ere taken from 
Woosley & Weaver ( [l995| ) and Maed er (|l992j ). SNIa data are 
from Nomoto, Thielmann fc Yok oi (1984). For intermediate 
mass stars Renzini & Voli (1981) data were used. The c ode 
also uses the SSPs of Weiss, Peletier & Matteucci (1995) to 
estimate the effects of non-solar abundance ratios on <Fe> 
and Mg2 predictions. Currently, our GCE code outputs esti- 
mates of 21 Lick indices for the population being modelled. 



4.1 Primordial Collapse. 

Primordial collapse represents the isolated collapse of a pri- 
mordial gas cloud. The cloud commences SF when the local 
density rises above some critical level. This mechanism has 
been proposed for the formation of elliptical galaxies. Dy- 
namical models of elliptical galaxy formation by primordial 



is < 1 G yr. The Chemo-dynamical models of of Theis et 
al. (1992) also show the main burst of star formation to be 



complete in <1 Gyr. Thus, these models predict a SFH com- 
prising a single burst of SF with a large majority of stars 
formed within the first 1 Gyr. This naturally leads to the 
prediction of [Mg/Fe] > solar, observed in the centres of el- 
liptical galaxies, as there is insufficient time for SNIa to con- 
tribute large quantities of Fe to the ISM before the bulk of 
SF is complete. These models also successfully predict many 
photometric properties such as the de Vaucouleurs r 1//4 pro- 
file and the presence of colour (metallicity) gradients. How- 
ever, several authors (Worthey et al. 1996; Greggio 1997; 
Vazdekis et al. 1997; Sansom & Proctor 1998) have shown 
that primordial collapse models fail to reproduce the strong 
central, metallicity sensitive, line-strengths found in ellipti- 
cal galaxies, under the assumption of a constant, Salpeter 
IMF. This is the result of the high number of low metal- 
licity stars produced in this scenario. We find this to be 
true also for the centres of spiral bulges which exhibit sim- 
ilar metallicities as the centres of elliptical galaxies (fig. ^). 
To demonstrate this point we have calculated indices for a 
range of primordial collapse models. Parameters within the 
models were set to maximise predicted central metallicities 
while keeping models realistic. Thus, as in-fall of primordial 
gas into the modelled region would dilute the metals con- 
tent of the star forming ISM, we assumed enriched inflow. 
The gas flow rate per Gyr was set to 10 times the initial 
gas content of the region. As the inflow is enriched, this 
large value ensures that the number of low metallicity stars 
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Figure 6. H/3 vs Fe4668 for SO galaxies (Fisher et al. 1996) and the 4 spiral bulges presented here. (N.B. Gorgas et al. (1997) did not 
include Fe4668.) Dotted lines indicate Worthcy's SSP values for metallicities; -0.25 0.0, +0.25, +0.5 and ages of 1.5 to 17 Gyr. Solid and 
dashed lines are as in fig. H (see section W). 



produced is kept small. Significantly higher values of inflow 
rate had negligible effects on the achieved line-strengths. To 
be consistent with dynamical models of primordial collapse 
the period of inflow was stopped after 0.4 Gyr while the SF 
was allowed to continue to the present day, albeit at an ever 
decreasing rate as gas is consumed. This generates an ~ 1 
Gyr burst of SF, which is also c onsist ent with the chemo- 
dynamical models of Theis et al. ( 1992 ). Continuation of SF 
to the present day also maximises the predicted strengths 
of metallicity dependent lines. To demonstrate the range of 
metallicities such models can achieve, star formation effi- 
ciency (C) during the inflow period was varied between 0.2 
and 4.0 Gyr -1 . This range of SF efficiencies s pans the val- 
ues given by Fritz-v. Alvensleben & Gerhard (1996) for Sb, 
Sa and E galaxies up to values for ultra-luminous star-burst 
galaxies. After the initial 0.4 Gyr of inflow, SF was allowed 
to continue in all models with an efficiency C=0.2. The ages 
of the models varied between 1.5 and 17 Gyr. The results of 
our GCE models of primordial collapse for <Fe>, Mg2, H/3 
and Fe4668 are shown as solid lines in figs. ^ to |(| Compari- 
son of indices for E, SO and spiral bulges (fig. |^) to the range 
of values achieved by primordial collapse models shows the 
inability of these models to reproduce the observed central 
line-strengths for both <Fe> and Mg2 features. This result 
is confirmed in figs | and §. Introduction to models of a bi- 
ased IMF at early epochs can resolve the under-abundance 



problem but can not simultaneously explain the high H/3 
values observed in some elliptical galaxies and at least one 
of our spiral bulges (fig. |B|). We might also expect to see an 
Mg2 excess with respect to solar in the biased IMF scenario. 
There is no evidence for this excess in our bulge data. 



4.2 Merger Models. 

While direct evidence for primordial collapse is not observed, 
examples of ong oing galaxy mergers are. Indeed, Schweizer 
& Seitzer ( 1993) suggest that, based on the number of galax- 



ies exhibiting on-going mergers or the fine structure indica- 
tive of a recent event, ~ 50% of field ellipticals have under- 
gone a merger event in the last 7 Gyr. This would seem to 
be supported by the fact that ~ 50% of the elliptical and SO 
sample shown in fig. [| show H/3 values that lie above the 8 
Gyr SSP line. In merger models, galaxy formation proceeds 
by coalescence of fragments that have undergone SF prior to 
merger. Such models can achieve high central line-strengths 
as the period of SF prior to merger pre-enriche s the ISM 
from which th e bulk of the population is formed ( sansom & 
Proctor 1998[) . It should be noted that both primordial col- 
lapse and merger models of galaxy formation predict short 
bursts of SF d ue to the extremely rapid inflow of gas i n 
both scenarios ( Theis et al. 1992 ; Barnes fc Hernquist 1996 ). 



© RAS, MNRAS 000, §-|l4| 



12 R. Proctor et al. 



<Fe> index (A) 



Fe4668 Index (A) 



Mg 2 Index (mag) 



H fl Index ( 





i 1 1 1 1 1 1 1 1 1 1 




0.5 1 
LOG(R) (arcsec) 



0.5 1 

LOG(R) (arcsec) 




1.5 



0.25 



0.2 



0.15 




0.5 1 1.5 

LOG(R) (arcsec) 



1.5 r 



0.5 



'.5 



0.5 



I I 




; ^ 




: ! 

m 



















0.5 1 1.5 

LOG(R) (arcsec) 



Figure 7. Index against log radius (radius in arcsec) for <Fe>, Fe4668, Mg2 and H/3 indices of the 4 galaxies of the Palomar data set. 
Open squares are disc contaminated side of the galaxies. All data points within the regions spanned by the regression lines are included 
in the fit. Note the flattening in metallicity sensitive indices in NGC 3623 which may be associated with the presence of a bar. 



Therefore, the excess Mg observed in elliptical galaxies is 
predicted by both of these models, without the need for the 
assumption of a biased IMF. The modelling of spiral bulges, 
using merger models, is constrained by their lack of Mg2 ex- 
cess. As the time-scale of gas inflow during an interaction is 
extremely rapid, large merger induced star-bursts will result 
in an Mg2 excess. As this is not seen in our spiral bulges, if 
bulge formation is dominated by accretions/mergers, the SF 
bursts must have been small and numerous. This is consis- 
tent with the presence of the disc which is highly unstable 
to me rgers of more than a f ew percent of the total galaxy 
mass (roth & Ostriker 1992). While we can construct single 
burst merger models that reproduce the <Fe>, H/3 and Mg2 
indices observed in the bulges, detailed modelling of merger 
scenarios has not been carried out as the large parameter 
space implied by multiple short bursts makes these models 
highly degenerate. In the limiting case, where interactions 
are so frequent that the induced SF is effectively continu- 
ous, we expect merger models converge to extended inflow 
models. 



4.3 Models with Extended Inflow. 

Extended inflow models represent long term inflow of gas 
into the region being modelled (e.g. Samland et al. 1997). 
The long duration of the inflow is consistent with formation 
of the bulge both by in-fall of gas from the disc (perhaps by 
the action of bars) as well as by the accretion of gas from the 
halo or captured satellite galaxies over an extended period. 



The models assume constant SF efficiency (C) with rela- 
tively modest gas inflow over a period of 2 to 17 Gyr and 
model spectral indices at a variety of points in the history 
of the population. With reasonable choices for C (0.4 - 4.0 
Gyr" 1 ) and inflow rate (10 s - 10 6 M Gyr -1 ) the observed 
<Fe>, Mg2 and H/3 in bulges can simultaneously be repro- 
duced to within errors. An inflow rate of 10 6 Mq Gyr -1 
corresponds to an increase in mass in the volume modelled 
equal to the initial mass every Gyr. The range of SF effi- 
ciencies excludes low values (0<C<0.4) as, even with con- 
tinuous SF for 17 Gyr, there is insufficient star formation to 
produce the metals required to achieve the high central line- 
strengths of NGC 3623 and NGC 4565. Formation episodes 
of 2 Gyr result in an Mg2 enhancement not reflected in 
the data (using the SSP models of Weiss et al. 1995). Figs. 
^ to ^| show the range of indices achieved by a selcetion of 
extended inflow models that reproduce the indices of the 
spiral bulges in both the Mg2 vs <Fe> and <Fe> vs H/3 
planes to within the errors. However, it should be noted 
that the models are highly degenerate with respect to in- 
flow rate and duration, SF efficiency, and time of SF onset 
(galaxy age). Consequently, the models are only shown to 
illustrate their ability to achieve the indices. The range of 
SF efficiencies and inflow durations that can reproduce the 
<Fe>, Mg2 and H/3 indices in our 4 bulges ar e con sistent 
with the findings of Samland, Hensler & Theis (1997) using 



a che mo-dy namical model of disc galaxy formation. Blum 
et al. (1996) showed that red giants in the Galactic bulge 



possess approximately solar metallicities and suggest that 
multiple epochs of SF have occurred in the centre of the 
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bulge in the last 7 - 100 Myr. Consequently, extended inflow 
models of spiral bulge formation are, so far, entirely consis- 
tent with both observations of Lick indices, the stars in our 
own bulge and chemo-dynamical modelling. 



5 CONCLUSIONS. 

Our study of the bulges of 4 spiral galaxies using the 
Lick system of spectral indices has shown that central line- 
strengths are high. <Fe> indices are similar to those found 
in the centres of elliptical galaxies. However, a difference in 
Mg2 index at a given <Fe> exists between elliptical galaxies 
and our sample of spiral bulges. Spiral bulges lie within the 
region of the Mg2 vs <Fe> pla ne oc cupied by the solar abun- 
dance ratio SSPs of Worthey (1994), while ellipticals exhibit 
enhanced Mg2 (the known [Mg/Fe] excess). This difference 
between the two object types reflects differences in their star 
formation histories. Our data are consistent with both cor- 
relations of Mg2 with central velocity dispersion and Mg2 
gradient with central velocity dispersion observed in ellipti- 
cal galaxies. 

Using our GCE code, we have shown that the central 
line-strengths in spiral bulges cannot be achieved with pri- 
mordial collapse models of spheroid formation with or with- 
out the assumption of constant, Salpeter IMF. Indeed, the 
inferred solar [Mg/Fe] argues against a biased IMF in spiral 
bulges. Models of bulge formation with gas inflow and star 
formation extended over 2-17 Gyr can achieve the observed 
central line-strengths in all of our sample. These models are 
consistent with t he che mo-dynamical modelling of Samland, 
Hensler & Theis (1997) as well as obser vations of individual 



stars in the bulg e of our own Galaxy (Minniti et al. 1992 



Blum et al. 1996 ). It has been shown that at least one bulge 



(NGC 5689) must contain a population of relatively young 
stars (^ 5 Gyr) again consistent with extended inflow mod- 
els of bulge formation. 

NGC 3623 shows the presence of kinematic substruc- 
ture. A dip in velocity dispersion is observed in the centre 
of this galaxy coincident with a sharp peak in metallicity 
sensitive indices. This structure suggests the presence of a 
bar, or perhaps disc, at the centre of this galaxy. Analysis of 
a much larger sample of galaxies, including a repeat obser- 
vation of NGC 3623, is being carried out with observations 
from the WHT, to test the findings of this paper. These 
data will be presented in a future paper and will include 
other age sensitive indices, such as the newly calibrated H5 
and H7 indices, which are less affected by emission than the 
H/3 index. Increased accuracy of age determination, coupled 
with the larger number of metallicity sensitive indices, will 
permit tighter constraint of our models of bulge formation 
and their star formation histories. 
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